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Исследуются возмущения ионосферной плазмы, создаваемые газовой струей реактивного двигателя космического аппарата, способные 
оказывать существенное воздействие на распространение и рассеяние радиоволн. Анализируются механизмы образования возмущений 
ионосферной плазмы космическим аппаратом с работающим двигателем. Излагается стохастическая модель процесса диффузии ионов 
ионосферной плазмы внутрь газовой струи реактивного двигателя космического аппарата. Приводятся результаты расчетов образования 
области высокой концентрации заряженных частиц ионосферной плазмы в окрестности границы газовой струи реактивного двигателя косми- 
ческого аппарата. На основе проведенного анализа предлагаются модели плазменных возмущений, создаваемых газовой струей реактивного 
двигателя космического аппарата в ионосфере. 


ІЙ?5Щ§эвание влияния возмущений ионосферной 
плазмы, создаваемых запускаемыми космическими 
аппаратами (КА), на распространение радиоволн отно- 
сится к числу важнейших проблем физики ионосферы 
и распространения радиоволн. Решение указанной 
проблемы необходимо для изучения процессов, про- 
исходящих в верхних слоях атмосферы при наличии 
таких возмущений, и выработки рекомендаций для 
устойчивой работы систем радиосвязи, радиолокации 
и навигации. 

Цель данной работы состоит в исследовании про- 
цессов создания ионосферных возмущений факелом 
КА, способных оказывать существенное воздействие 
на распространение и рассеяние радиоволн. 

Начиная с первых запусков мощных ракет, спо- 
собных выводить на орбиту искусственные спутники 
Земли, проблеме взаимодействия газовых струй 
реактивных двигателей КА с верхними слоями ат- 
мосферы уделяется большое внимание. Истечение 
газовых струй в затопленное пространство и вакуум, 
а также их взаимодействие со спутным потоком 
набегающего газа рассматривалось в ряде работ 
[1-5]. 

Несмотря на большое количество публикаций, 
посвященных перечисленным вопросам, в настоящее 
время нельзя сделать вывод о том, что проблема взаи- 
модействия газовых струй реактивных двигателей КА 
с верхней атмосферой решена. Трудности построения 


адекватных моделей взаимодействия обусловлены как 
сложностью расчета всех процессов, происходящих в 
факеле ракеты [5], так и необходимостью учета влияния 
неоднородности атмосферы [4], а также геомагнитного 
поля на формирование структуры возмущения концен- 
трации заряженных частиц [6]. 

Широкий спектр проблем, связанных с взаимо- 
действием истекающей газовой струи КА с верхней 
атмосферой, не позволяет рассмотреть их все. Поэтому 
ниже проводится анализ основных типов ионосферных 
возмущений, вызываемых газовой струей КА, способ- 
ных оказывать существенное воздействие на распро- 
странение радиоволн. Основное внимание уделяется 
разработке модели воздействия КА на ионосферную 
плазму в области ближней зоны до расстояния, на ко- 
тором происходит выравнивание давлений газа струи 
и атмосферы. 

1. Механизмы возмущения ионосферной плазмы газовой 
струей реактивного двигателя космического аппарата 

Истечение газов из сопла реактивного двигателя на 
высотах 100... 400 км от поверхности Земли характери- 
зуется большим перепадом давлений на срезе сопла (р ц ) 
и в ионосфере (р 0 ). Согласно [5, 7], давление на срезе 
сопла для жидкостных реактивных двигателей состав- 
ляет величину р ц = 0,02... 0, 5 кг/см 2 , а для двигателей на 
твердом топливе р ц = 1 ...3,5 кг/см 2 

Полагая р ц = 1 кг/см 2 и используя данные о давлении 
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М 


в ионосфере, приведенные в [8], можно определить 
пределы изменения степени нерасчетности п на высотах 

100... 400 км... я ._ , ... 

ѵ = \/“ о =3-10 б ...3-10 ш . 

— "исел Маха на срезе сопла 

. ,. /и =2. ..5 г , ,и. 

( ’ - направленная и тепловая 

скорости исм "его газа) для жидкостных реактивных 

двигателей и ' _ для двигателей на твердом топливе 
[5, 7, 9]. ЧисЛд л Іаха набегающего спутного потока при 
скорм — У . п “7, ° '"Э мо) изменяться в преде- 
лах 0_ 0 и ‘ - ■” ’ где 11 - тепловая скорость 

частиц в ионосфере. 

Согласно [5, 7], степень иор т '™ ,т '" ""за на срезе 
сопла изменяется в диапазоне 


10 ... 10 в 
ферной плазме степеш ті '' ,,ті уации изменяется от 


= 10- 


на высоте 100 км до 


= 10 


на высоте 400 км. Длина 


свободного пру 1 """ "еитральных частиц в и<у'~~'*'зре 
возрастает от = * см на высоте 100 км до = ^ км 
на высоте 400 км ПО-12]. 

На высотах ^ км магнитное поле Земли начина- 

ет оказывать существенное влияние на движение ионов 
плазмы, поскольку длина свободного пробега частиц 

гг е ларморовского радиуса ионов 

. >р №). = 2...4 


Для строгого решения задачи взаимодействия га- 
зовой струи реактивного двигателя КА с ионосферой 
необходимо совместное решение системы кинетических 
уравнений, записанных для всех сортов частиц взаимо- 
действующих газов факела и ионосферы. Трудности 
усугубляются присутствием в ионосфере магнитного 
поля, влияние которого ня пкижение ионов необходимо 
учитывать на высотах - ' км [6]. 

Поскольку в такой постановке решение задачи 
практически невозможно, то в ряде работ по взаи- 
модействию газовой струи реактивного двигателя 
с ионосферой полагается, что при описании взаи- 
модействия применимы уравнения гидродинамики, 
а ионы и электроны увлекаются нейтральными 
частицами (роль магнитного поля не учитывается) 
[5,7,13]. 

В работах [1-4] в указанном приближении проведе- 
но численное решение задачи о расширении сверхзвуко- 
вых струй в вакуум [1, 2], в затопленное пространство 
[1, 3] и в спутный сверхзвуковой поток набегающего 
газа [1, 4]. 

Вследствие вытеснения ионосферной плазмы 
слабоионизованным газом струи (степень ионизации 
которого, как правило, меньше ионосферной на высотах 
- 100 км) на границе струи образуется область высокой 
концентрации заряженных частиц ионосферной плазмы 
(ОВК), а внутри струи - область низкой концентрации 
заряженных частиц (ОНК) [6]. 

С другой стороны, газовую струю можно рассмат- 
ривать как некоторое пористое тело, движущееся в ио- 
носфере со скоростью КА. При движении газовой струи 
происходит как отражение частиц ионосферной плазмы 
от струи, так и проникновение их внутрь. Отраженные от 
струи частицы, сталкиваясь с равновесными частицами 
ионосферной плазмы, возмущают их концентрацию и 
температуру. Аналогичный случай рассмотрен в работе 


[14], где источником отраженных молекул является по- 
верхность твердого тела. В результате возмущение обра- 
зует ударную волну, которая способна распространяться 
в ионосферной плазме, возмущая ее концентрацию. 

Ионы ионосферной плазмы, проникают внутрь 
струи, сталкиваются с частицами газовой струи, теряя 
при этом энергию и импульс, вследствие чего происхо- 
дит их торможение и образование ОВК. Столкновения 
между частицами газовой струи и ионосферы приводят 
к повышению температуры газа в области смешения. 
Часть энергии при столкновениях может переходить во 
внутреннюю энергию молекул газа. В результате воз- 
можна дополнительная термоионизация газа в области 
смешения. Подобная задача рассмотрена в работе [15] 
для энергичных молекул, отраженных от поверхности 
движущегося КА. В следующем разделе остановимся 
подробнее на результатах численного моделирования 
процесса образования ОВК и ОНК. 

2. Моделирование процесса диффузии ионов ионосферной 

плазмы внутрь факела космического аппарата 

При сверхзвуковом движении слабоионизованной 
газовой струи реактивного двигателя КА в ионосфере на 
высотах = 200 км ионы ионосферной плазмы могут 
проникать внутрь газовой струи, отдавая при столк- 
новениях свою энергию и импульс частицам струи, и 
увязать в ней, образуя область повышенной концентра- 
ции заряженных частиц. Решение задачи проводится в 
приближении химически не реагирующих газов. Плазма 
предполагается квазинейтральной, так что возмущения 
концентрации электронов и ионов совпадают. 

Модель газовой струи выбирается, исходя из закона 
сохранения числа частип расширяющегося слоя газа 
N Г Г 

(сііѵ ( ) = 0, где N и - концентрация и скорость 

частиц газа). В цилиндрической системе координат, 
связанной с источником (КА), где ось направлена в 
противоположную сторону относительно направления 
движения КА, для макропараметров газа струи полу- 
чаются выражения 

Рі 7\ = о, Р / ? >і ё ѳ 0 ,^, 
т м /\= с о7^р/^ѳ 0 , 

т = сош!, 1 ” ? = СОП5І, ^ = ? ? р 7, 

С С 

где Ѳ 0 - угол полураствг^ — " р " - координаты 
начала и конца струи, о — " / і > - радиус среза 

сопла, ]Ч 0 - концентрация нейтральных частиц струи 
на срезе сопла, Т, Ѵ 2 и Ѵ р - температура и компоненты 
скорости истекающего газа. 

Задача о диффузии ионов разреженной плазму у) 
движущуюся со сверхзвуковой скоростью струю, задат 
ваемую моделью ( 1-3), решается методом Монте-Карло 
[16]. Для однородного газа в разделе 1.2 монографШ 
[6] излагается схема построения дерева возмущения, 
вызываемого пробной частицей в поле частиц, распре- 
деленных по максвелловскому закону. Ниже указанная 
схема применяется для построения ветви возмущения, 
вызываемого ионом, движущимся в неоднородном 
газовом образовании. 

При моделировании на ПЭВМ координатное про- 
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странство делится на ячейки, ог-рниченные поверх- 
ностями Ѳ. = соші и плоскостями = сопй. Предпола- 
гается, что отношение значений концентрации частиц 
•у _ N /Ы ’ ячейках постоянно и близко к 1, 

1 — ' При этом координаты границ ячеек 
растут, как э г; _ у 72 : ( 2 ) по закону геометрической 
прогрессии ~с " с, ' Г~РУ Я ограничивается по оси 
плоскостью - “ ’ " - выбирается из условий: 

превышения газового давления в струе к давлению в 
ионосфере и малого отношеюу —--бега иона в 

к-ой ячейке к размеру ячейки ' 1 ' ‘ 

По углу Ѳ границы ячеек выбираются равномерно, 
начиная с некоторого Ѳ . и до Ѳ . В оперативной памяти 
ПЭВМ выделяется соответствующий массив ячеек для 
вычисляемой величины возмущения концентрации ио- 
нов Д1Ч. Перед началом работы программы элементам 
этого массива присваиваются нулевые значения. 

Потоку ионов ионосферной плазмы, набегающему 
на струю, при численном моделировании процесса 
диффузии ставится в соответствие источник ионов, рас- 
положенный на поверхности струи. Выбор начальных 
координат и скорости иона проводится стандартным 
способом моделирования произвольного распределения 
Р(у) случайной величины у, изложенным в Приложении 
1 монографии [ 6 ]. 

Вероятность нахождения иона на поверхности 
струи в интервале координаты определяется 
выражениеіѵ л; <; \ы; 


П/5\6С = . 


>46 Л) 2 


поперечной оси 


-2*6 


і соотношением 


Л = Ф(= , )+ ехр 6 ° ” 2 )/(Ѵ* ? ѵ ^ 

? ѵ = ? 0 5ІпѲ 0 , 

? 

0 - скорость движения КА, скорости и и V 

п и„ , 

ражены в единицах тепловой скорости ионов ’ 

ф 6)=7=і сх р(- Т2 Г (6) 

0 - интеграл вероятности.^) 

При вычислениях с целью сокращения времени 
счета выражения ( 6 - 8 ) табулируются. Выбор времени 
и места столкновения иона с нейтральной частипей 
факела проводится в системе координат, связанной с ло- 
ш 7 ?ѵ ЬІМ участком газа, перемещающимся со скоростью 
* * . При столкновении выполняются законы сохране- 
ния массы, импульса, энергии и числа входящих и вы- 
ходящих частиц. Предполагается, что взаимодействие 
между частицами происходит по закону жестких сфер. 
Используются выраженіи, применявшиеся в разделе 1 .2 
из [ 6 ] для нахождения времени между столкновениями 
и координат следующего столкновения. 

Скорости частиц су"С""’іи т„, т, и скоростями 

і_) п и і 

перед столкновением 05 •’ после столкновения по 


ЗАКОНУ У\ і -) _і_ I м -) ііѵіігі I і / I 

в;=-‘^+"|о 0 -б,|ц/(ц + і ), 

^ = Ч| | - У, -Ц||ц/(Ѵі+Ц 

Ц = М / М * „ „ 

где г 1 і / 2 » _ случайный вектор, равномерно 

распределенный на поверхности сферы единичного 
радиуса, при моделировании которого используются 
выражения (П1.8-П1.14) из [ 6 ]. 

При слежении за к-ьв ^ іном проводится накопление 
времени его пребывания в фой ячейке координатного 
пространства. Наблюдение за ионом прекращается при 
выходе его за границы струи. После построения М вет- 
вей значение относительного возмущения концентрации 
ионов в ’’ — -\е координатного пространства (при 

ч ~ .. , 

условии ) для начальной максвелловской фун- 

кции распределения ионов по скопостям определяется 
выражением [ 6 ] 


О 0 +К 


ми [17] 


Из (4) и (П1.5) из [ 6 ] следует выражение для коор- 
динаты рожд рииа 

^0+6 Л )>/6 

где х 2 - случайное число, равномерно распределенное 
в интервале [ 0 , 1 ]. 

Азимутальный угол ф, отсчитываемый в плоскости 


где х^ - случайное число, равномерно распределенное 
в интервале [ 0 , 1 ]. (4) 

Вероятность-ппобной частице (иону) иметь скорость 
и в интервале а поверхности струи, движущейся 
со скоростью КА 0 относительно ионосферной плаз- 
мы, впредь- --- 8 г-ч5 - 
(и) и = х (и+% у (и) и, 

V. 

где - норму" ~ поверхности струи, с - константа 
нормировки, ' ' - функция распределения частиц 
по скоростям. В результате подстановки (5) в (П1.5) из 
[ 6 ] в локальной системе координат с осью ’ направ- 
ленной по нормали к поверхнос ( . получаются 

трансцендентные птпштшчп.пп ’ уравнения 

\=(1 + Ф(иОУ2, 

^=(і + ф(о,)Ѵ2, 

У - ехр (-и? У (Ѵт? , )+ Ф (и? )} 

3 = (Г76 ’ (5) 


где °‘ -равновесная концентрация ионов ионосферной 
плазмы, 8 - площадь боковой поверхности газовой струи, 
\Ѵ. - объем ]-ой ячейки координатного пространства, 
М - число ветвей возмущения, по которым проводится 
усреднение. 

Затем весь процесс повторяется и при увеличении 
числа ветвей М результаты, полученные при различных 
реализациях, стремятся к одному значению. После 
усреднения по всем реализациям получается искомое 
среднее значение величины относительного возмуще- 
ния концентрации ионов. Следует отметить, что вели- 
чина относительной ошибк” т'^р~тений уменьшается 
с ростом числа ветвей как 1 ' Объем статистики 

при расчетах в данной работе выбирался таким, чтобы 
обеспечить относительной погрешнос® 

вычислений - 1 • 
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Результаты проведенных расчетов представлены на 

рис. 1-4 для сл<^"~”~ значений параметров: скорость 

движения КА 0 - ’ ско — с™ т, и в системе 
координат, связанной с КА, и - и,ѵ ( и ' - тепловая 
скорость нейтральных частиц газа струи); отношение 
концентрации нейтральных частиц в соседних по г 
ячейках § п = 1 , 1 ; отношение длины свободного пробега 
нейтрал <; _ р , ~ аза стр \ ,Ѵ _ .г сопла к ве- 

не » _ '7 1 ё°о равно ' " = г° ’ ’) и 

1 ГГ 6 Г ь /Ч _ 1 ^АА ’ 

ш (рис. 4); для этих же случаев “ > 0 и 

4000 соответственно; отношение іц — ' к массе 


41КД7 соответственно; отношение іц к массе 

ней™"""™' """тицы газ?-- _ ' _І (рис. 1, 

4Ѵ 7 " (рис. 2), 7 ѵ - Гпис. На 


(рис. 3). На 


верхней половине рисунков изо^~й";ц" ’ изолинии, а на 

А / ' , 

нижнеи- величина возмущения ’ 0 в зависимости 
от координаты р, отсчитываемой от оси струи. 

Анализ полученных результатов показывает, что 
максимум относительного возмущения концентрации 
ионов располагается в области, прилега’"^' АГ ~~~-це 
струи, и составляет величину порядка л ' °' = ' 

Эффект накопления ионов объясняется торможе- 
нием ионов в газовой струе за счет их столкновений с 
нейтральными частицами газа струи и увеличением их 
времени жизни в ячейках координатного пространства. 
Это, согласно (9), и приводит к образованию области 
уплотнения. 

С уменьшением о^-ул — ия массы иона к массе 
нейтральной частицы ' ' (рис. 1-3) размеры воз- 
мущенной области увеличиваются, что объясняется тем, 
что легкие ионы в течение более длительного времени 
сохраняют энергию при столкновениях с массивными 
нейтральными частицами и дальше проникают внутрь 
струи. 

Уменьшение отношения длины свободного пробе- 
га частиц газовой струи 1 на срезе сопла к величине 


5 _ Р 

' ! ® 0 соответствует увеличению ксу — ™" 

ции і. ", V л?'"™' газа струи. Та ѵ "Г 

= 10 . і о 6 ■ =10 18 .3 

для ' 0 N |9 и получается см 3 

(рис. 1-3), '=10 см -з (р ИС 4) соответственно. Уве- 

личение плотности газовой струи приводит (рис. 1, 4) 
к возрастанию числа столкновений ионов с нейтраль- 
ными частицами струи и к росту величины максимума 
относительного возмущения концентрации ионов. Как 
следует из сравнения результатов, представленных 
на рис. 1 и 4, увеличение расхода топлива на порядок 

у^г ИТ к росту размеров возмущенной области в 

= ѵЮ раз. 

3. Модели возмущений ионосферной плазмы, создаваемых газовой 

струей реактивного двигателя космического аппарата 

На основе проведенного анализа механизмов возму- 
щения ионосферной плазмы газовой струей реактивного 
двигателя КА ниже выбираются модели плазменного 
образования, создаваемого КА в ионосфере [6]. Рас- 
сматривается область ближней зоны до расстояния, на 
котором происходит выравнивание давлений химически 
не взаимодействующих газа истекающей струи и атмос- 
ферного газа. 

Первая модель работает на высотах -- ,<)() км, 
рис. 5 (I). В этом случае, вследствие большой длины 
свободного пробега нейтральных частиц ионосферной 
плазмы — 1 км, возмущение концентрации заряженных 
частиц, вызываемое нейтральный и, отра- 
женными от газовой струи, мало: Л ' 0 < В то же 

время механизм увязания ионов ионосферной плазмы 
в газовой струе реактивного двигателя КА работает на 
этих высотах и возмущение концентрации заряженных 
частиц ионосферной плазмы представляет собой (как 
это показали результаты предшествующих расчетов) 
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Естественные науки 



дІМ /Ы 

Рис. 3. Относительное возмущение концентрации ионов ' ' 



Д М /N 

Рис. 4. Относительное возмущение концентрации ионов ' ' 



г 


осесимметричное квазиконическое плазменное обра- 
зование. 

Относительное значение возмущения концентрации 
заряженных частиц ионосферной плазмы в началь- 
ном участке струи может достигать величины = 5- 
Коэффициент отражения радиоволн от поверхности 
плазменного образования определяется в основном 
концентрацией электронов и частотой их столкновений 
с нейтральными молекулами. Оценки показывают, что 
в среднем в обз 4 ия частота электронных 

столкновений ѵ ‘ - ' В эт <|> < лч - ттз ф е- 

'К N “К Ф ,лн с частотой _ V / « ( 

_ о 1 + ^ » - плазменная частота ионосферы 

на высоте движения КА) коэфф ф отражения от 

возмущения близок к 1, так как 4 ' [18, 19]. 

Известно, что при работе двигателя КА струя исте- 
кающих газов турбулизована. Источником турбулент- 
ностей могут служить как неравновесные процессы в 
камере сгорания, сопровождаемые флуктуациями давле- 
ния и температуры, так и неустойчивости, возникающие 
при взаимодействии нескольких струй газа, истекающих 
в пространство с малым давлением. 

Вследствие этого, область уплотнения концентрации 
заряженных частиц ионосферной плазмы, создаваемая 
факелом КА, будет также турбулизована. Наличию 
турбулентных вихрей в плазменном образовании можно 
поставить в соответствие элементарные рассеиватели 
[20, 21], перемещающиеся вдоль образующей пла^' '~н- 
ного образования со скоростью истечения струи " ‘ 

Вторая модель справедлива до высот ионосферы 

- ' 30 км, где длина свободного пробега газе 

сравнима с ларморовским радиусом ионов = Р 1 " - рис. 
5 (II). В этом случае влиянием магнитного поля на дви- 
жение ионов можно пренебречь и в начальном участке 
струи возмущение представляет собой совокупность 
внутреннего и внешнего плазменных образований. 
Внутреннее плазменное образование (ОВК), как и в 
первой модели, обусловлено работой механизма увя- 
зания ионов. Образование ударной волны частицами, 
отраженными от факела КА, которое возможно на этих 
высотах, вследствие малой длины свободного пробега 
частиц газа, приводит к появлению внешнего плазмен- 
ного образования. 

Третья моделію зуется на высотах ионосферы 

130. . .300 км, где - Р и необходимо учитывать вли- 
яние магнитного поля на формирование возмущения 
заряженных частиц ионосферной плазмы, вызываемого 
ударной волной, генерируемой факелом КА, рис. 5 
(III). 



] М о.- \ 


ш 


Рис. 5. Модели плазменного образования 


Заключение 

1 . Проведен анализ механизмов создания возмуще- 
ний ионосферной плазмы газовой струей реактивного 
двигателя КА. 

2. Исследован механизм образования области вы- 
сокой концентрации заряженных частиц ионосферной 
плазмы газовой струей реактивного двигателя КА. 
Результаты расчетов показывают, что максимум отно- 
сительного возмущения концентрации ионов располага- 
ется в области, прилегающей к границе струи. Эффект 
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накопления ионов объясняется торможением ионов в 
газовой струе за счет их столкновений с нейтральными 
частицами плотного газа струи. 

3. На основе проведенного анализа механизмов воз- 
мущения ионосферной плазмы газовой струей выбраны 
модели плазменных образований, создаваемых факелом 
КА в верхней атмосфере. 
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